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Abstract 

We study here a natural situation when constraint programming can be entirely reduced 
to rule-based programming. To this end we explain first how one can compute on con- 
straint satisfaction problems using rules represented by simple first-order formulas. Then 
we consider constraint satisfaction problems that are based on predefined, explicitly given 
constraints. To solve them we first derive rules from these explicitly given constraints and 
limit the computation process to a repeated application of these rules, combined with 
labeling. 

We consider here two types of rules. The first type, that we call equality rules, leads to a 
new notion of local consistency, called rule consistency that turns out to be weaker than arc 



consistency fo r constraints of arbitrary arity (called hyper-arc consistency in (Marriott & 
^tuckey, 1998 )). For Boolean constraints rule consistency coincides with the closure under 



the well-known propagation rules for Boolean constraints. The second type of rules, that 
we call membership rules, yields a rule-based characterization of arc consistency. 

To show feasibility of this rule-based approach to constraint programming we show how 



both types of rules can be automatically generated, as CHR rules of (Friihwirth, 1995). This 
yields an implementation of this approach to programming by means of constraint logic 
programming. 

We illustrate the usefulness of this approach to constraint programming by discussing 
various examples, including Boolean constraints, two typical examples of many valued 
logics, constraints dealing with Waltz's language for describing polyhedral scenes, and 
Allen's qualitative approach to temporal logic. 



Note. A preliminary version of this article appeared as ( Apt fc Monfroy, 199£ ). In this 
version we also present a framework for computing with rules on constraint satisfaction 
problems and discuss in detail the results of various experiments. 



1 Introduction 
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1 . 1 Background 

This paper is concerned with two styles of programming: constraint programming 
and rule-based programming. 

In constraint programming the programming process is limited to a generation 
of constraints and a solution of the so obtained constraint satisfaction problems 
(CSP's) by general or domain dependent methods. 

In rule-based programming the programming process consists of a repeated ap- 
plication of rules. A theoretical basis for this programming paradigm consists of 



so-called production rules that were introduced in the seventies, see, e.g., (Luger 
fc Stubblefield, 199^ ) [pages 171-186], though the idea goes back to the works of 
A. Thue and of E. Post in first half of twentieth century. The production rules 
are condition-action pairs, where the condition part is used to determine whether 
the rule is applicable and the action part defines the action to be taken. The most 
known programming language built around this programming paradigm was 0PS5 



of ( iForgy, 198l| ) 



Recently, there has been a revival of interest in rule-based programming in the 
context of constraint programming. The earliest example is the CHR language of 



( Friihwirth, 1995 ) that is a part of the ECL'PS'^ system. (For a more recent and 
more complete overview of CHR see ( [Friihwirth, 1998| ).) The CHR rules extend the 
syntax of constraint logic programming by allowing two atoms in the conclusion and 
employing guards. These rules are predominantly used to write constraint solvers. 

Another example of a programming language in which rules play an important 
role is ELAN. It offers a logical environment for specifying and prototyping deduc- 
tion systems by means of conditional rewrite rules controlled by strategies. ELAN 
is used to support the design of various rule-based algorithms such as constraints 
solvers, decision procedures, theorem provers, and algorithms expressed in logic 
programming languages, and to provide a modular framework for studying their 



combinations. A general overview of ELAN can be found in (Borovansky et 



1998 ), whereas ( Kirchner fc Ringeissen, 199 j ) and ( Castro, 1998 ) (to which we 
shall return in Section |ll|) describe applications of ELAN to constraint program- 
ming and constraint solving. 

Also, in the hybrid functional and object-oriented language programming lan- 
guage CLAIRE of (^ aseau &: Laburthe, 1996| ) rules are present. CLAIRE was de- 
signed to apply constraint programming techniques for operations research prob- 
lems. The rule-based programming is supported by means of production rules that 
can be naturally used to express constraint propagation. 

It is useful to mention here that also logic programming and constraint logic 
programming are rule-based formalisms. However, these formalisms use rules dif- 
ferently than rule-based programming described above. This distinction is usually 
captured by referring to forward chaining and backward chaining. In rule-based 
programming, as discussed above, forward chaining is used while in logic program- 
ming and constraint logic programming backward chaining is employed. Intuitively, 
forward chaining aims at a simplification of the considered problem and it main- 
tains equivalence, while backward chaining models reasoning by cases, where each 
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case is implicitly represented by a different rule. Both forms of chaining can be 
combined and in fact such a combination is realized in the CHR language, in which 
the CHR rules model forward chaining while the usual Prolog rules model backward 
chaining. 

1.2 Overview of Our Approach 

The traditional way of solving CSP's consists of combining constraint propaga- 
tion techniques with search. Constraint propagation aims at reducing a CSP to an 
equivalent one but simpler. In case of finite domains the most basic approach to 
search consists of labeling, a repeated enumeration of the domains of the successive 
variables. 

The aim of this paper is to show that constraint programming can be entirely 
rendered by means of rule-based programming. To this end we provide a framework 
in which one computes on CSP's by means of rules represented by simple first- 
order formulas. In this approach the constraint propagation is achieved by repeated 
application of the rules while search is limited to labeling. This yields a framework 
for constraint programming more related to logic than the usual one based on 
algorithms achieving local consistency. 

The rules we shall consider are implications built out of simple atomic formulas. 
In our study we focus on two types of rules. The first type, that we call equality 
rules, are of the form 

Xl = . . .,Xn = Sn ^t 

where a;i, . . ., Xn, y are variables and si, . . ., s„, t are elements of the respective vari- 
able domains. The computational interpretation of such a rule is: 

if for i € [l..n] the domain of the variable Xi equals the singleton {si}, then remove the 
element t from the domain of y. 

The second type of rules, that we call membership rules, are of the form 

Xl G Si, . . .,Xn € Sn ^ y ^ t 

where 

variables and Si,. . .,Sn are subsets of the respective variable 

domains, 

• y is a variable and t is an element of its domain. 
The computational interpretation of such a rule is: 

if for i e [l..n] the domain of the variable x, is included in the set Si, then remove the 
element t from the domain of y. 

To illustrate the use of these rules we study CSP's that are built out of predefined, 
explicitly given finite constraints. Such CSP's often arise in practice. Examples in- 
clude Boolean constraints, constraints dealing with Waltz's language for describing 
polyhedral scenes, Allen's temporal logic, and constraints in any multi- valued logic. 

To solve such CSP's we explore the structure of these explicitly given constraints 
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first. This information is expressed in terms of valid equality and membership rules. 
The computation process for a CSP built out of these constraints is consists of two 
phases: a generation of the rules from the explicitly given constraints and a repeated 
application of these rules, combined with labeling. 

To characterize the effect of the generated equality and membership rules we use 
the notions of local consistency. These notions approximate in a loose sense the 
notion of "global consistency", see, e.g., ( Tsang, 1993 ). We show that the first type 
of rules leads to a local consistency notion that turns out to be weaker than arc 
consistency for constraints of arbitrary arity. We call it rule consistency. 

When the original domains are all unary or binary, rule consistency coincides with 
arc consistency. When additionally the predefined constraints are the truth tables 
of the Boolean connectives, these rules are similar to the well-known Boolean prop- 
agation rules, (see, e.g., ( Friihwirth, 1998 ) [page 113]). As a side effect, this shows 
that the Boolean propagation rules characterize arc consistency. Rule consistency 
is thus a generalization of the Boolean propagation to non-binary domains. 

We also show that the membership rules lead to a notion of local consistency 
that coincides with arc consistency. This yields a rule-based implementation of arc 
consistency. 

To show feasibility of this rule-based approach to constraint programming, we 
automatically generate both types of rules, for an explicitly given finite constraint, 
as rules in the CHR language. When combined with a labeling procedure such CHR 
programs constitute automatically derived decision procedures for the considered 
CSP's, expressed on the constraint programming language level. In particular, we 
automatically generate the algorithms that enforce rule consistency and arc consis- 
tency. 

The availability of the algorithms that enforce rule consistency and arc con- 
sistency on the constraint programming language level further contributes to the 
automation of the programming process within the constraint programming frame- 
work. In fact, in the case of such CSP's built out of predefined, explicitly given finite 
constraints, the user does not need to write one's own CHR rules for the considered 
constraints and can simply adopt all or some of the rules that are automatically 
generated. In the final example of the paper we also show how using the equality 
rules and the membership rules, we can implement more powerful notions of local 
consistency. 

Alternatively, the generated equality rules and membership rules could be fed 
into any of the generic Chaotic Iteration algorithms of (Apt, 1999a) and made 
available in such systems as the ILOG solver. This would yield rule consistency and 
an alternative implementation of arc consistency. 

The algorithms that for an explicitly given finite constraint generate the appro- 
priate rules that characterize rule consistency and arc consistency have (unavoid- 
ably) a running time that is exponential in the number of constraint variables and 
consequently are in general impractical. 

To test the usefulness of these algorithms for small finite domains we implemented 
them in ECL^PS*^ and successfully used them on several examples including the ones 
mentioned above. The fact that we could handle these examples shows that this 
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approach is of practical value and in particular can be used to automatically derive 
practical decision procedures for constraint satisfaction problems defined over small 
finite domains. Also it shows the usefulness of the CHR language for an automatic 
generation of constraint solvers and of decision procedures. 



1.3 Organization of the Paper 

The rest of the paper is organized as follows. In the next section we clarify the 
syntax of the rules and explain how one can compute with them. In Section 3 we 
illustrate the use of these computations by means of an example. In Section 1 we 
prove that the outcomes of the computations we are interested in are unique. In 
Section ^ we introduce some semantic aspects of the rules and in Section ^ we 
formalize the concept of a CSP is built out of predefined constraints. Next, in 
Section |^ we introduce the notion of rule consistency and discuss an algorithm that 
can be used to generate the minimal set of rules that characterize this notion of local 
consistency. Then, in Section || we compare rule consistency with arc consistency. 
In Section ^ we study membership rules and discuss an algorithm analogous to the 
one of Section |^. This entails a notion of local consistency that turns out to be 
equivalent to arc consistency. 

In Section [l^ we discuss the implementation of both algorithms. They gener- 
ate from an explicit representation of a finite constraint a set of CHR rules that 
characterize respectively rule consistency and arc consistency. We also illustrate 
the usefulness of these implementations by means of several examples. Finally, in 
Section ^ we discuss other works in which a link was made between constraint 
programming and rule-based programming and in Section |l^ we assess the merits 
of our approach. In the appendix we summarize the tests carried out by means of 
our implementation of both algorithms. 



2 Computing with Rules 

In what follows we introduce specific type of rules and explain how one can compute 
with them on constraint satisfaction problems. First, we introduce constraints. 

Consider a sequence of variables X :~ xi, . . .,Xn where n > 0, with respective 
domains Di, . . .,Dn associated with them. So each variable Xi ranges over the do- 
main Di. By a constraint C on X we mean a subset of Di x . . . x D„. Given an 
element d :— di, . . dn of Di x . . . x I?„ and a subsequence Y :— Xi-^ , . . ., Xig. of X 
we denote by d\Y] the sequence di-^, . . di,. In particular, for a variable Xi from X, 
d[xi] denotes di. 

Next, we define the rules we are interested in. 

Definition 2.1 

• Let X be a variable, a an element and S a set. By an atomic formula we mean 
one of the following formulas: x = a, x^a, x€iS. 

• By a rule we mean an expression of the form Ai, . . Am — > -Bi, . . ., -B„, where 
each Aj and is an atomic formula. □ 
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In what follows a rule will be always associated with some constraint. Then every 
atomic formula x = a or x ^ a (respectively, x ^ S) will be such that a belongs to 
the domain of x (respectively, 5 is a subset of the domain of x) . 

Subsequently, wc explain how to compute using the rules in presence of con- 
straints. First, we limit our considerations to the rules of the form Ai, . . A^ — » a; ^ 
a. We need to explain how to turn the disequality formula into an action. This is 
done by identifying the disequality x ^ a with the assignment ■= — {a}, 
where is the current domain of x. In other words, we interpret a; 7^ a as an 
action of removing the value a from the current domain of the variable x. 

This leads us to the definition of an application of such a rule. We need some 
semantic notions first. 

Definition 2.2 

Consider a constraint C on a sequence of variables X, a variable x of X, and a 
tuple d & C. 

• Given an atomic formula A involving x we define the relation \=d A as follows: 

— |=d a; = a iff d[x] = a, 

— |=d a; 7^ a iff d[x] ^ a, 

— hd .X e 5 iff d[x] e S. 

• Given a sequence of atomic formulas A := Ai,. . .,Am we define \=d A iff 
\=d Ai for all i G [l..m]. □ 

Definition 2.3 

Consider a constraint C on a finite sequence of variables X and a rule of the form 
A Xi a involving only variables from X. 

Suppose that for all d G C wc have \=d A. Let C" be the constraint obtained 
from C by removing the element a from the domain Z)j of the variable Xi and by 
removing from C all tuples d such that d[xi] = a. Then we call the constraint C 
the result of applying the rule A — > a;^ 7^ a to C. 

If a G Di, then we say that this is a relevant application of the rule A — > ^ a 
to C. If C" coincides with C, we say that this application of the rule A — > aij ^ a 
to C maintains equivalence. □ 

So the application of the rule A — > aij 7^ a to a constraint C on the sequence 
a;i, . . ., a;„ of variables with respective domains Z>i, . . ., £)„ results in the constraint 
C" on the variables a;i, . . ., a;„ with respective domains Di, . . ., -Dj-i, .D^, -Di+ij • • D„ 
where 

• Di = Di- {a}, 

• C = C n (£>i X . . . X A-i xD'iX Di+i, . . ., D„). 

We say then that the constraint C is restricted to the domains Di, . . Z)j_i, £)■, 
-Dj+l, . . ., Dn- 

Now that wc defined the result of a single application of a rule we proceed to 
define computations. To this end we first introduce constraint satisfaction problems. 
By a constraint satisfaction problem, in short CSP, we mean a finite sequence of 
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variables X with respective domains V, together with a finite set C of constraints, 
each on a subsequence of X. We write it as {C ; xi G Di, . . .,Xn 6 where 
X := xi, . . .,Xn and V := Di, . . ., Dn- 

By a solution to (C ; xi S -Di, . . ., a;n G -D„) we mean an element d G Di x . . . x I?„ 
such that for each constraint C € C on a sequence of variables X we have d[X] G C. 
We call a CSP consistent if it has a solution. Two CSP's with the same sequence 
of variables are called equivalent if they have the same set of solutions. 

We now modify the definition of an application of a rule to a constraint to an 
application of a rule to a CSP. To this end we attach each rule to a constraint 
to which it is supposed to be applied. Even though the constraints change during 
the computations wc consider, it will be always clear from the context to which 
constraint a given rule is attached. 

Definition 2.4 

Consider a CSP V and a rule Xi ^ a attached to a constraint C of V. Suppose 
that for all d G C we have \=d A. Define a CSP V on the same variables as V as 
follows: 

• the domain of Xi in V equals Di — {a}, where Di is the domain of Xi in V, 

• the domains of other variables in V' are the same as in V, 

• the constraints of V' are obtained by restricting the constraints of V to the 
new domains. 

We say then that the CSP V is the result of applying the rule A ^ Xi ^ a to V. 
U a £ Di, then we say that this is a relevant application of the rule A — > ajj 7^ a to 

v. □ 

Finally, we introduce the crucial notion of a computation. 
Definition 2.5 

Consider a set of rules 7?. of the form A^ x ^ a and an initial CSP V. By a 
computation by means of TZ starting at V we mean a maximal sequence of CSP's 
Vi, . . .,Vi, . . . such that each Vj+i is the result of a relevant application of a rule 

bomUtor-j. □ 

Note that when the set of rules TZ is finite or when all domains in P are finite, all 
computations starting at V are finite. The reason is that in each of these two cases 
the number of elements mentioned in the conclusions of the rules in TZ is finite. 
But each element can be removed from a domain only once and we insist that in 
computations each rule application is relevant, from which the claim follows. If a 
computation is finite, then no application of a rule from the considered set of rules 
TZ to the final CSP is relevant, i.e., this final CSP is closed under the rules in TZ. 

The computations are a means to reduce the domains of the variables while 
preserving the equivalence of the considered CSP. The computations here considered 
are in general insufficient for solving a CSP and in the case of CSP's with finite 
domains they have to be combined with labeling. Labeling can be modeled in the 
above rule-based framework by introducing a rule that splits a given CSP into two, 
the union of which is equivalent to the given CSP. The addition of such a rule to 
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the considered framework leads to no conceptual difficulties and is omitted. On the 
other hand, various forms of labeling strategies, like the one in which variable with 
the smallest domain is chosen first, cannot be captured on this level. 

The above string of definitions allowed us to define computations in which the 
actions are limited to the applications of rules of the form A — > a; 7^ a acting on 
CSP's. 

By limiting our attention to such type of rules we do not lose any expressiveness. 
Indeed, consider first a rule of the form A ^ a; = a. To compute with it we interpret 
the equality a; = a as the assignment n {a}, where is the current 

domain of x. Then each rule A — > x = a is equivalent to the conjunction of the 
rules of the form A x ^ b with b £ D — {a}, where D is the original domain of 

X. 

Next, consider a rule of the form A ^ a; G S*. To compute with it we interpret 
the atomic formula x € S as the assignment :— n S, where is the current 
domain of x. Then each rule A ^ a; G 5 is equivalent to the conjunction of the 
rules of the form A — > x 7^ 6 with 6 G -D — S*, where D is the original domain of x. 

Finally, each rule of the form A ^ Bi, . . Bm is equivalent to the conjunction of 
the rules of the form A^ Bi for i G [l..m]. 

Note that the rules of the form xi = oi, . . ., a;„ = an ^ y ^ b are more expressive 
than so-called dependency rules of database systems (see, e.g., ( UUman, 1988| )). 
These are rules of the form a;i = oi, . . ., a:„ = a„ y = b. We just explained how 
to model them by means of rules of the form xi = oi, . . ., a:„ — an ^ y ^ b. 

However, modeling in the other direction is not possible, as can be seen by taking 
the variables x,y, each with the domain {0,1,2}, and the constraint C on x,y 
represented by the following table: 





y 





1 








2 


2 



Then the rule x 
rules. 



^ y 7^ 2 is not equivalent to a conjunction of the dependency 



3 An Example 

We now show how we can use the rules for computing by means of an example 
kindly provided to us by Victor Marek. We solve here a simple logic puzzle from 
( [Fleming, 2doo| ). 

Below, given a set of variables a;i, . . ., a;„, each with the domain D we denote the 
following set of rules: 



{xi = a ^ Xj ^ a \ i,j G [l..n],i / j, a G D} 
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by alLdifferent{xi, . . ., These rules formalize the requirement that the variables 
all different. The puzzle is as follows. 

To stave off boredom on a rainy Saturday afternoon, Ms. Rojas invented a game for 
Denise and her two other children to play. Each child selected a different household object 
(no two of which were in the same room) to describe to the others, who would try to guess 
the item and its location in the house. Can you match each child with the item he or she 
selected to describe, as well as the room of the house (one is the living room) in which each 
is located? 

Here arc the clues provided: 

1. The three children are Byron, the child who selected the book, and the one 
whose item is in the den; 

2. The rug is in the dining room; 

3. Felicia selected the picture frame. 

To solve this puzzle we use nine variables, 

• childi, child2, childs, to denote the three children, Byron, Denise and Felicia, 

• roomB,roomD,roomF, to denote the rooms of, respectively, Byron, Denise 
and Felicia, 

• item B, item D, item F, to denote the objects selected by, respectively, Byron, 
Denise and Felicia. 

We postulate that the domain of childi is {Byron}, of child2 is {Denise} and 
of child^ is {Felicia}. Next, we assume that each rocmii variable has the set 
{den, dining, living} as its domain and each itemi variable has the set {book, frame, rug} 
as its domain. 

The initial set up of the story is formalized by the following rules: 

• all_different{roomB, roomu, roomp), 

• alLdifJerent{itemB, itemo, itemp)- 

This yields 36 rules but we shall group the rules with the same premise, so we 
shall actually have 18 rules. The rules we shall need below will be, from the first 
set: 

(rl) rocrniB = dining roomo 7^ dining, roomp 7^ dining, 
(r2) roomB = living roomo 7^ living, roomp 7^ living, 
(r3) roomo = living roomB 7^ living, roomp ^ living, 

and from the second set: 

(r4) itemB = rug itemn ^ rug, itemp ^ rug, 

(r5) itemp = frame — > itemB 7^ frame, itemo 7^ frame. 

Next, the first clue is formalized by means of eight rules out of which the only 
ones of relevance below will be 

(cl.l) itemB 7^ book, 

(cl.2) itemo = book — > roomo 7^ den. 
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The second clue is formalized by means of six rules out of which the only one of 
relevance below will be 

(c2) iteniB = rug roairiB 7^ den, rooruB 7^ living. 

Finally, the third clue is formalized by means of two rules: 

(c3.1) — > itemp 7^ rug, 
(c4.1) — > itemp 7^ book. 

So in total we have 34 rules but we shall use below only the ten rules made 
explicit. The initial CSP has nine variables as introduced above and one single 
"universal" constraint that consists of the Cartesian product of all the variable 
domains. The computation consists of twelve steps and proceeds as follows. 

1. Using the rule (cS.l) the domain of itemp is limited to {book, frame}. 

2. Using the rule (c3.2) the domain ofitemp is further limited to {frame}. Thus 
itemp = frame is established. 

3. Using the rule (r5) and the fact itemp = frame just established the domain 
of iterriB is limited to {book, rug}. 

4. Using the rule (cl.l) the domain of iterriB is further limited to {rug}. Thus 
iteniB = rug is established. 

5. Using the rule (r5) and the fact itemp = frame established in step 2 the 
domain of itemo is limited to {book, rug}. 

6. Using the rule (r4) and the conclusion items = rug of step 4 the domain of 
itemo is further limited to {book}. Thus itemo = book is established. 

7. Using the rule (c2) and the fact items = rug established in step 4 the domain 
of rooms is limited to {dining, living}. 

8. Again using the rule (c2) and the fact items = rug established in step 4 the 
domain of rooms is further limited to {dining}. Thus rooms = dining is 
established. 

9. Using the fact rooms = dining established in step 8 and the rule (rl) the 
domain of roomD is limited to {den, living}. 

10. Using the fact itemp, = book established in step 6 and the rule (cl.2), the 
domain of roomo is further limited to {living}. Thus roomo = living is 
established. 

11. Using the fact rooms = dining established in step 9 and the rule (r2) the 
domain of roomp is limited to {den, dining}. 

12. Using the fact roomo = living established in step 10 and the rule (r3) the do- 
main of roomp is further limited to {den}. Thus roomp = den is established. 

At this stage one can check that the resulting CSP with all singleton domains is 
closed under all 34 rules. This yields the solution to the puzzle represented by the 
following table: 
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child 


room 


item 


Byron 


dining 


rug 


Denise 


living 


book 


Felicia 


den 


frame 



4 Outcomes of Computations 

A natural question arises whether the outcome of computations using a finite set 
of rules is unique. The answer is positive. To prove it we need a lemma concerning 
iterations of inflationary and monotonic functions. 

Definition 4-1 

Consider a partial ordering {D, □ ) with the least element _L and a finite set of 
functions F -.^ {fi,.. ., fk} on D. 

• By an iteration of F we mean an infinite sequence of values do, di, . . . defined 
inductively by 

do := ±, 

dj := ^.(dj-i), 

where each ij is an element of [l..k]. 

• We say that an increasing sequence do ^ di ^ d2 . . . of elements from D 
eventually stabilizes at d if for some j > we have di = d for i > j. 

• A function / on Z? is called inflationary if x Q fix) for all x. 

• A function / on D is called monotonic if a; C y implies f{x) C f{y) for all 
x,y. □ 



We now need the following lemma from ( Apt, 19991: ) . 



Lemma 4-2 

Consider a partial ordering (D, C ) with the least element _L and a finite set of 
monotonic functions F on D. Suppose that an iteration of F eventually stabilizes 
at a common fixpoint d of the functions from F. Then d is the least common fixed 
point of the functions from F. 

It follows that all iterations of F that eventually stabilize at a common fixpoint 
stabilize at the same element. We now prove the desired result. 



Theorem 4-3 

Fix an initial CSP V. Consider a finite set TZ of rules of the form A 
all computations by means of TZ starting at V yield the same CSP. 



X ^ a. Then 



Proof 
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We already noted in Section ^ that all such computations arc finite. Suppose that 
V := {C ; xi G Di, . . .,Xn G Dn). We consider now the following partial ordering 
C ). The elements of D are the sequences (£'1, . . ., En) such that Ei C Di for i 6 
[l..n], ordered componentwise w.r.t. the reversed subset ordering 3. So (Z?i, . . ., £)„) 
is the least element _L in this ordering and 

{Eu...,En) C iSE.DF, for i £ [l..n]. 

We replace in each rule each premise atom Xi = a hy Xi £ {a} and Xi ^ a hy 
Xi G Di — {a}. Since for all d G Di x . . . x Z?„ we have \=d 2;^ = a iff \=d Xi G {a} 
and [^d Xi ^ a iff \^d G Di — {a}, it follows that the applications of the original 
and of the resulting rules coincide. This allows us to confine our attention to the 
rule each premise of which is of the form z G S. 

Consider now a membership rule zi G 5*1, . . ., Zm G Sm y =/= o, associated with 
a constraint C from C defined on a set of variables Y. We interpret this rule as a 
function on the just defined set D as follows. 

First, denote by C the extension "by padding" of C to all the variables xi, . . ., a;„, 
i.e. C C Di X . . . X Dn and d e C* iff d[Y] e C. Next, given a constraint E and its 
variable z denote the set {d[z] \ d G E} by E[z]. Finally, assume for simplicity that 
y is Xn- 

The function / that corresponds to the rule zi e 5i, . . ., Zm G Sm y ^ t 
defined as follows: 

r (-El, . . ., En ~ {a}) if (C n X . . . X En))\z^\ C 
/(El,... ,£;„):= <^ foriG[l..m], 
[ . . ., i?„) otherwise. 

Denote the set of so defined functions by F . By definition each function / G F is 
infiationary and monotonic w.r.t. the componentwise reversed subset ordering 3. 

Now, there is a one-one correspondence between the common fixpoints of the 
functions from F at which the iterations of F eventually stabilize and the outcomes 
of the computations by means of TZ starting at P. In this correspondence a common 
flxpoint (El, . . ., En) is related to the CSP (C ; xi G i?i,...,a;„ G En) closed 
under the rules of TZ, where C are the constraints from C restricted to the domains 



El, . . .,En- The conclusion now follows by Lemma 4.2. □ 



5 Semantic Aspects of Rules 

We now introduce a number of semantic notions concerning rules. 

Definition 5.1 
Consider a constraint C. 

• We say that the rule A ^ B is valid for C if for all tuples d £ C 

A imphes 1=^ B. 

• We say that the constraint C is closed under the rule A — > B if 

{^d A for all tuples d G C) implies {\=d B for all tuples d G C). 
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• We say that the rule A — > B is feasible for C if for some tuple d G C we have 
N A. 

• We say that the rule A ^ B for the constraint C extends the rule A' ^ B if 
A contains all variables of A' and for all tuples d € C 

\=d A implies \=d A'. 

• Given a set of rules TZ, we call a rule minimal in TZ if it is feasible and it does 
not properly extend a valid rule in TZ. □ 

To illustrate them consider the following example. 

Example 5.2 

Take as a constraint the ternary relation that represents the conjunction and{x, y, z). 
It can be viewed as the following table: 



a; I y I ^; 



In other words, we assume that each of the variables x, y, z has the domain {0, 1} 
and view and{x, y, z) as the constraint on x, y, z that consists of the above four 
triples. 

It is easy to see that the rule a; = 0-^2;7^1is valid for and{x, y, z). Further, 
the rule x = Q,y = 1 ^ z ^ 1 extends the rule x = — > 2; 7^ 1 and is also valid for 

and{x, y, z). However, out of these two rules only x — 0—>-z^lis minimal. 

Finally, both rules arc feasible, while the rules x = 0,z = l^yj^O and x = 

0, 2=1— >2/7^1 arc not feasible. □ 

Note that the definition of an application of a rule is so designed that a link with 
semantics is kept in the following sense: if a rule r is valid for a constraint C, then 
C is closed under r. Rules that are not feasible are trivially valid. Note also that a 
rule that extends a valid rule is valid, as well. So validity extends "upwards" . 

Note the use of the condition "A contains all variables of A'" in the definition 
of the relation "the rule A ^ B extends the rule A' — > B" . Without it we would 
have the following paradoxical situation. Consider the variables x, y, z, all on the 
domains {0, 1}, and the constraint C on x, y, z defined by C := {(0, 0, 0), (1, 0, 1)}. 
Then the rules x = ^ z ^ 1 and y = ^ z ^ 1 are both valid for C and for all 
d G C we have that \=d a; = implies \=d j/ = 0. So without the mentioned condition 
we would have that the rule x = ^ z ^ 1 extends the rule y = ^ z 1, which 
would imply that the first rule is not minimal. 

In the sequel the following observation will be useful. 



Note 5.3 
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Consider two finite and non-empty constraints C and E such that C C E and a set 
of rules TZ. Then C is closed under all valid rules from TZ for E iff it is closed under 
all minimal valid rules in TZ for E. 

Proof 

Suppose that C is closed under all minimal valid rules in TZ for E. Take a rule r 
from TZ that is valid for E. 

Case 1. r is feasible for E. 

Then, because E is finite, r extends some minimal valid rule r' in TZ for E. But 
C is closed under r', so it is closed under r, as well. 

Case 2. r is not feasible for E. 

Then r is not feasible for C either since C E. Consequently, since C is non- 
empty, C is closed under r. □ 

In what follows we confine our attention to computations involving two types of 
rules: 

• equality rules: these are rules of the form xi = si, . . ., Xm = Sm y a; we 
abbreviate them to X — s y ^ a, where X = xi, . . Xm and s — si, . . .Sm, 

• membership rules: these are rules of the form xi £ Si, . . Xm G Sm —^y^a] 
we abbreviate them to X £ S ^ y ^ a, where X — xi,...,Xm and S = 

Si, . . .Sm- 



By specializing in the last clause of Definition 5.1, defining a minimal rule, the 
set TZ of rules to the set of equality rules and to the set of membership rules we 
obtain the notions of a minimal equality rule and of a minimal membership rule. 
For equality and membership rules the following straightforward characterization 
of the "extends" relation will be of use. 

Note 5.4 

i. An equality rule xi = si,...,Xm = Sm^y ^ a extends an equality rule 
zi = ti,...,Zn = tn ^ y a iff zi = ti, . . ., Zn = tn is a subsequence of 

Xi — Si, . . ., Xm — -^m- 

ii. A membership rule xi £ 5*1, . . ., Xm G Sm ~^ U ^ o, extends a membership rule 
zi e Ti, . . ., z„ & Tn y ^ a iff zi, . . ., z„ is a subsequence of xi, . . ., Xm and 
for each i G [\..ri\ we have S'^(i) C Ti, where Zi equals x^^i^iy □ 

Given a CSP with finite domains we would like to solve it by considering com- 
putations starting at it. But where do we get the rules from? Note that given a 
constraint C and a rule r that is valid for C, the constraint C is trivially closed 
under r. Consequently, an application of r to C is not relevant, i.e., it does not 
affect C. So to obtain some change we need to use rules that are not valid for the 
initial constraints. This brings us to the notion of a CSP based on another one. 

6 CSP's Built out of Predefined Constraints 

In the introduction we informally referred to the notion of a CSP "being built out 
of predefined, explicitly given constraints." Let us make now this concept formal. 
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We need two auxiliary notions first, where in preparation for the next definition we 
already consider constraints together with the domains over which they are defined. 

Definition 6.1 

• Given a constraint C C Di x . . . x _D„ and a permutation tt of [l..n\ we denote 
by the relation defined by 

(ai, . . ., a„) e C iff (a^(i), . . ., a^(„)) e C 

and call it a permutation of C . 

• Given two constraints C C £)i x . . . x _D„ and £^ C Z?'j^ x . . . x 1?^^ we say that 
C is based on E if 

— AC D[ for i e [1..7i], 

— C = EC^{DlX ...X Dn). □ 

So the notion of "being based on" involves the domains of both constraints. If G is 
based on E, then C is the restriction of E to the domains over which C is defined. 

Definition 6.2 

We assume that the "predefined constraints" are presented as a given in advance 
CSP BASE. Suppose that each constraint C of a CSP V is based on a permutation 
of a constraint of BASE. We say then that V is based on BASE. □ 

In the above definition the use of permutations of constraints allows us to abstract 
from the variable ordering used in BASE. The following example illustrates this 
notion. 

Example 6.3 

Consider the well-known full adder circuit. It is defined by the following formula: 

add{ii,i2,i3, 01:02) = 
xor{ii,i2,Xi),and{ii,i2, ai), xor{xi,i3, 02), and(i3, xi, 02), or{ai, 02, oi), 

where and, xor and or are defined in the expected way. We can view the original 
constraints as the following CSP: 

BOOC := (and{x, y, z), xor(x, y, z), or{x, y,z) ; x G {0, 1}, y G {0, 1}, z G {0, 1}}. 

BOOC should be viewed just as an "inventory" of the predefined constraints and 
not as a CSP to be solved. Now, any query concerning the full adder can be viewed 
as a CSP based on BOOC For example, in Section |l^ we shall consider the query 
add{l, x, y, z, 0). It corresponds to the following CSP based on BOOC: 

{ xor{ii,i2,Xi), and{ii,i2,ai), xor{xi,i3, 02), and(i3, xi, 02), or{ai, a2, oi) ; 
n G {l},i2 G {0,1}, Z3 G {O,l},0i G {0,1}, 02 G {0},ai G {0, 1}, 02 G {0,1}, 
XI G {0,1}). 

□ 

In what follows we consider computations that start with a CSP based on some 
CSP BASE. In these computations we wish to maintain equivalence between the 
successive CSP's. To this end the following simple observation is crucial. 
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Note 6.4 

Consider two constraints C and E such that C is based on E. Let A ^ a; 7^ a be a 
rule vahd for E. Then the apphcation of A ^ .t ^ a to C maintains equivalence. 

Proof 

Assume that the rule A ^ x 7^ a can be applied to C, i.e., that for all d G C we have 
\=d A. Suppose now that the rule A — > cc 7^ a to C does not maintain equivalence. 
Then for some d ^ C we have d[x] = a. C is based on E, so d ^ E. By the validity 
of the rule for E we get d[x] ^ a. This yields a contradiction. □ 

This observation provides us with a way of maintaining equivalence during a 
computation: it suffices to use at each step a rule that is valid for a permutation 
of a constraint of BASS. Such a rule is then attached (i.e., applied) to the 
constraint based on . This is what we shall do in the sequel. Depending on the 
type of rules used we obtain in this way different notions of local consistency. 



7 Rule Consistency 

In this section we consider a CSP V based on some finite CSP BASE and study 
computations that use exclusively equality rules. The rules are obtained from the 
constraints of BASE; each of them is valid for a permutation of a constraint 



C of BASE and is applied to the constraint of V based on . By Note 6.4 the 
successive CSP's are all equivalent to the initial CSP V. The computation ends 
when a CSP is obtained that is closed under the rules used. This brings us to a 
natural notion of local consistency expressed in terms of equality rules. 

Definition 7.1 

Consider a CSP V based on a CSP BASE. Let C be a constraint of V. For some 
constraint /(C) of BASE and a permutation tt, C is based on f{Cy. 

• We call the constraint C rule consistent (w.r.t. BASE) if it is closed under 
all equality rules that are valid for J{CY . 

• We call the CSP V rule consistent (w.r.t. BASE) if all its constraints are rule 
consistent. □ 

In what follows we drop the reference to BASE if it is clear from the context. 

Example 7.2 

Take as the base CSP 

BASE := {and{x, y,z) ; xe {0, 1}, y G {0, 1}, z G {0, 1}) 

and consider the following four CSP's based on it: 



1. {and{x,y,z) 

2. {and{x,y,z) 

3. {and{x,y,z) 

4. {and(x,y,z) 



x&{0},y<EDy,ze{0}), 
xG{l},yG7^„zG{0,l}), 
X G {0,l},y G Dy,ze {1}}, 
x&{0},yeDy,ze{0,l}). 
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where Dy is a subset of {0, 1}. We noted in Example 5.2 that the equality rule 
X — Q z ^ 1 is valid for and{x,y, z). In the first three CSP's its only constraint 
is closed under this rule, while in the fourth one it is not closed since 1 is present 
in the domain of z whereas the domain of x equals {0}. So the fourth CSP is not 
rule consistent. One can show that the first two CSP's are rule consistent, while 
the third one is not rule consistent since it is not closed under the valid equality 
rule z = 1 ^ X ^ 0. 



When trying to generate all valid equality rules Note 5.3 allows us to confine our 
attention to the minimal valid equality rules. We now introduce an algorithm that 
given a finite constraint generates the set of all minimal valid equality rules for it. 
We collect the generated rules in a list. We denote below the empty list by empty 
and the result of insertion of an element r into a list L by insert (r, L). 

By an assignment to a sequence of variables X we mean here an element s from 
the Cartesian product of the domains of variables of X such that for some d G C 
we have d[X] = s. Intuitively, if we represent the constraint C as a table with rows 
corresponding to the elements (tuples) of C and the columns corresponding to the 
variables of C, then an assignment to X is a tuple of elements that appears in some 
row in the columns that correspond to the variables of X. This algorithm has the 
following form, where we assume that the considered constraint C is defined on a 
sequence of variables VAR of cardinality n. 

Equality Rules Generation algorithm 

L : = empty ; 

FOR i:= TO n-1 DO 

FOR each subset X of VAR of cardinality i DO 
FOR each assignment s to J DO 
FOR each y in VAR-X DO 

FOR each element d from the domain of y DO 
r := X = s y ^ i; 

IF r is valid for C 

and it does not extend an element of L 
THEN inserter, L) 

END 

END 

END 

END 

END 

END 

The test that one equality rule does not extend another can be easily implemented 
by means of Note 5.4.i. 



The following result establishes correctness of this algorithm. 
Theorem 7.3 

Given a constraint C the Equality Rules Generation algorithm produces in L 
the set of all minimal valid equality rules for C. 
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Proof 

First note that in the algorithm aU possible feasible equality rules are considered 
and in the list L only the valid equality rules are retained. Additionally, a valid 
equality rule is retained only if it does not extend a rule already present in L. 
Finally, the equality rules are considered in the order according to which those 



that use less variables are considered first. By virtue of Note 5.4.i this implies that 
if a rule r2 extends a rule ri, then ri is considered first. As a consequence precisely 
all minimal valid equality rules are retained in L. □ 

The above algorithm is completely straightforward and consequently inefficient. 
It is easy to see that given a constraint defined over n variables, 0(n-2" ■d'^) rules are 
considered in it, where d is the size of the largest variable domain. This shows that 
in practice this algorithm is impractical for large domains and for constraints with 
many variables. By representing the rules explicitly one could improve the running 
time of this algorithm, trading time for space. Then the test that one rule does not 
extend another could be eliminated from the algorithm by representing explicitly 
the partial ordering defined by the relation "equality rule ri extends equality rule 
r2" . Each time an equality rule that is valid for C would be found, all the rules that 
extend it would be now disregarded in further considerations. This would reduce 
the number of rules considered and improve the average running time. However, 
it is difficult to quantify the gain obtained and in the worst case still all the rules 
would have to be considered. 

In Section |l^ and the appendix we present some empirical results showing when 
the Equality Rules Generation becomes infeasible. 



8 Relating Rule Consistency to Arc Consistency 

To clarify the status of rule consistency we compare it now to the notion of arc 



consistency. This notion was introduced in ( Mackworth, 1977a ) for binary relations 
and was extended to arbitrary relations in ( Mohr fc Masini, 1988 ). Let us recall the 
definition. 

Definition 8.1 

• We call a constraint C on a sequence of variables X arc consistent if for every 
variable a; in A" and an element a in its domain there exists d G C such that 
a = d[x]. That is, each element in each domain participates in a solution to 
C. 

• We call a CSP arc consistent if all its constraints are arc consistent. □ 

The following result relates for constraints of arbitrary arity arc consistency to 
rule consistency. 

Theorem 8.2 

Consider a CSP V based on a CSP BASE. If V is arc consistent, then it is rule 
consistent w.r.t. BASE. 



Proof 
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Assume that V is arc consistent. Choose a constraint C of 'P and consider an 
equahty rule X — s ^ y ^ a that is vaUd for f{CY, where / and tt are as in 



Definition 7.1 



Suppose by contradiction that C is not cfosed under this rule. So for X := 
Xi, . . .,Xk and s := si, . . ., the domain of each variable xj in V equals {sj} and 
moreover a ^ D, where D is the domain of the variable y vnV. 

By the arc consistency of V there exists d G C such that d[y] — a. Because of 
the form of the domains of the variables in X, also d[X] — s holds. Additionally, 
because V is based on BASE, we have d G f{CY. But by assumption the equality 
rule X = s ^ y a \s valid for f{CY, so d[y] ^ a. A contradiction. □ 

The converse implication does not hold in general as the following example shows. 

Example 8.3 

Take as the base the following CSP 

BASE := (C ; x G {0, 1, 2}, y e {0, 1, 2}) 

where the constraint C ox\ x,y that equals the set {(0, 1), (1,0), (2,2)}. So C can 
be viewed as the following table: 





y 





1 


1 





2 


2 



Next, take for Di the set {0,1} and D2 the set {0,1,2}. Then the CSP (C n 
(-Di X i?2) ; X e y e D2), so ({(0, 1), (1, 0)} ; x G {0, 1}, y G {0, 1, 2}) is based 
on BASE but is not arc consistent since the value 2 in the domain of y does not 
participate in any solution. Yet, it is easy to show that the only constraint of this 
CSP is closed under all equality rules that are valid for C. □ 

We now show that if each domain has at most two elements, then the notions of 
arc consistency and rule consistency coincide. More precisely, the following result 
holds. 

Theorem 8.4 

Let BASE be a CSP each domain of which is unary or binary. Consider a CSP V 
based on BASE. Then V is arc consistent iff it is rule consistent w.r.t. BASE. 



Proof 

The ( ) implication is the contents of Theorem 3.2. 

To prove the reverse implication suppose that some constraint C of T' is not arc 
consistent. We prove that then C is not rule consistent. 

The constraint C is on some variables xi, . . ., a;„ with respective domains Di, . . ., Z)„ 
For some i G [l..n] some a € Di does not participate in any solution to C. 



20 



Krzysztof R. Apt and Eric Monfroy 



Let Di-^ , . . ., Dii, be the sequence of all domains in Di, . . ., -Di_i, i^i+i, • ■ •, £)„ that 
are singletons. Suppose that Di. :— {si- } for j G and let X :— Xi^, . . ., Xi^ and 

S . Si-^^ , . . ., Si^ . 

Consider now the equality rule X = s ^ Xi ^ a and take /(C)'^, where / and 
TT are as in Definition 7.f . For appropriate domains D'j^, . . of BASE we have 
/(C)- C i?; X . . . X D'^. 

Next, take some c? e ]{GY such that = s. We show that d € C. Since 

C = f{CY n (Di X ... X £)„) it sufhces to prove that d € Di x . . . x Z?„. For each 
variable Xj lying inside of X we have d[xj\ — Sj € Dj. In turn, for each variable 
Xj lying outside of X its domain Dj has two elements, so, by the assumption on 
BASE, Dj is the same as the corresponding domain D'j of f{CY and consequently 
d[xj] e Dj, since d & D[ x . . . x D'^. 

So indeed d ^ C and hence o?[a;i] 7^ a by the choice of a. This proves validity of 
the equality rule X — s ^ Xi a ioi f{CY- 

But C is not closed under this rule since a e Z^i, so C is not rule consistent. □ 



9 Membership Rule Consistency 

In this section we consider computations that use exclusively membership rules. In 
the previous section we saw that the notion of rule consistency is weaker than that 
of arc consistency for constraints of arbitrary arity. Here we show that by using the 
membership rules we obtain a notion of local consistency that coincides with arc 
consistency. 

First, let us clarify the notion of a membership rule by considering the following 
example. 

Example 9.1 

Consider a constraint on variables x, y, z, each with the domain {+, — , I, r}, that is 
defined by the following table: 



a; I 2/ I z 



+ 


+ 


+ 


I 


7' 






I 


r 


r 




I 



This constraint is the so-called fork junction in the language of ( Waltz, 1975 ) for 
describing polyhedral scenes. Note that the following three membership rules 

ri X e {+, -} z ^ I, 
r2 ■■— X e {+} ~> z ^ I, 
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and 

rs-.^xe {-},y e{l}'^z^l 

are all valid. The membership rules r2 and r^ extend ri while the membership 
rule ri extends neither r2 nor ri. Further, the membership rules ^2 and r^, are 
incomparable in the sense that none extends the other. □ 



Now, in analogy to Definition 7.1, we introduce the following notion. 



Definition 9.2 

Consider a CSP V is based on a CSP BASE. Let C be a constraint of V. For some 
constraint /(C) of BASS and a permutation tt, C is based on f{CY. 

• We call the constraint C membership rule consistent (w.r.t. BASE) if it is 
closed under all membership rules that are valid for f{CY. 

• We call a CSP membership rule consistent (w.r.t. BASE) if all its constraints 
are membership rule consistent. □ 

We now have the following result. 

Theorem 9.3 

Consider a CSP V based on a CSP BASE. Then V is arc consistent iff it is mem- 
bership rule consistent w.r.t. BASE. 

Proof 

( ^ ) This part of the proof is a simple modification of the proof of Theorem B.2 . 

Assume that V is arc consistent. Choose a constraint C ol V and consider a 
membership rule X £ S ^ y ^ a that is valid for f{C)'^, where / and tt are as in 
Definition 



6.2 



Suppose by contradiction that C is not closed under this rule. So for X 
xi, . . .,Xk and 5 := 5i, . . ., S'fc the domain of each variable Xj is included in Sj and 
moreover a G D, where D is the domain of the variable y. 

By the arc consistency of V there exists d G C such that d[y] — a. Because of 
the form of the domains of the variables in X, also d[xi] e Si for i G [l..k] holds. 
Additionally, because V is based on BASE we have d € f{C)^. But by assumption 
the rule X € S ^ y ^ a is valid for f{C)^, so d[y] ^ a. A contradiction. 



( 4= ) This part of the proof is a modification of the proof of Theorem 3.4. 

Suppose that some constraint C of T' is not arc consistent. We prove that then C 
is not membership rule consistent. The constraint C is on some variables xi, . . .,Xn 
with respective domains iDi,...,Z3„. For some i £ [l..n] some a G Di does not 
participate in any solution to C. 



Take f{C)'^, where / and tt are as in Definition 9.2. For appropriate domains 
£);,...,£>; of BASE we have /(C)^ C D[ x . . . x D'^. 

Let Di^ , . . ., Dii be the sequence of domains in _Di, . . ., -D^-i, . . ., Dn that 

are respectively different than D'l, . . ., D[_i, -D^+i, . . ., -D^j. Further, let A" := x^j^ , . . .,Xi 



and S := Di^,.. .,Dii,. 



Consider now the membership rule A e 5 ^ / a. Take some d G f{C)^ such 
that d[x,^] e A, for j G [1..^]. We show that d G C. Since C = f{C)-^ n (Di x 
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. . . X Dn) it suffices to prove that d G Di x . . . x £>„. For each variable Xj lying 
inside of X we have d[xj] G Dj. In turn, for each variable Xj lying outside of X 
its domain Dj is the same as the corresponding domain of f{CY in BASE and 
consequently d[xj] G Dj, since d G D'-^^x . . . x Z?^. 

So indeed d G C and hence (i[a;i] ^ a by the choice of a. This proves validity of 
the rule X £ S ^ Xi ^ a ior f{CY. But C is not closed under this membership 
rule since a G so C is not membership rule consistent. □ 



Example 8.3 shows that the notions of rule consistency and membership rule 
consistency do not coincide. To see this difference better let us reconsider the CSP 
discussed in this example. 

We noted there that this CSP is not arc consistent and that it is rule consistent. 
From the above theorem we know that this CSP is not membership rule consistent. 
In fact, consider the following membership rule: 

X e {0, l}^y^ 2. 

This membership rule is valid for the base constraint C but the restricted con- 
straint C n (Di X D2) is not closed under this rule. In conclusion, the membership 
rules are more powerful than the equality rules. 

As in Section |^ we now provide an algorithm that given a constraint generates 
the set of all minimal valid membership rules. We assume here that the considered 
constraint C is defined on a sequence of variables VAR oi cardinality n. 

Instead of assignments that are used in the Equality Rules Generation 
algorithm we now need a slightly different notion. To define it for each variable 
X from VAR recall that we denoted the set {d[x] \ d G C} by C[x]. By a weak 
assignment to a sequence of variables X := Xi,. . .,Xk we mean here a sequence 
Si, . . ., S'fc of subsets of, respectively, C[xi], . . ., C[xk] such that some d G C exists 
such that dlxi] G Si for each i e [l..k]. 

Intuitively, if we represent the constraint C as a table with rows corresponding 
to the elements of C and the columns corresponding to the variables of C and we 
view each column as a set of elements, then a weak assignment to X is a tuple 
of subsets of the columns that correspond to the variables of X that "shares" an 
assignment. 

In the algorithm below the weak assignments to a fixed sequence of variables are 
considered in decreasing order in the sense that if the weak assignments Si, . . .,Sk 
and Ui, . . ., Uk are such that for i € [l..k] we have Ui C Si, then Si, . . .,Sk is con- 
sidered first. 

Membership Rules Generation algorithm 

L : = empty ; 

FOR i:= TO n-1 DO 

FOR each subset X of VAR of cardinality i DO 

FDR each weak assignment S to X in decreasing order DO 
FOR each y in VAR-X DO 

FOR each element d from the domain of y DO 
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r := X e S y ^ d; 
IF r is valid for C 

and it does not extend an element of L 

THEN insert Tr, L) 

END 

END 

END 

END 

END 

END 

The test that one membership rule does not extend another can be implemented 



using Note 5.4.ii 



The following result establishes correctness of this algorithm. 
Theorem 9.4 

Given a constraint C the Membership Rules Generation algorithm produces 
in L the set of all minimal valid membership rules for C. 

Proof 



The proof is analogous to that of Theorem 7.3, We only need to check that the 



membership rules are considered in such an order that if a rule r2 extends a rule 



ri, then ri is considered first. This follows from directly from Note 5.4.ii. □ 



10 Applications 

In this section we discuss the implementation of the Equality Rules Genera- 
tion and Membership Rules Generation algorithms and discuss their use on 
selected domains. 



10.1 Constraint Handling Rules (CBR) 

In order to validate our approach we have realized in the Prolog platform ECL'PS'^ 
a prototype implementation of both the Rules Generation algorithm and the 
Membership Rules Generation algorithm. We made a compromise between 
memory usage and performance so that we could tackle some non-trivial problems 
(in terms of size of the domains of variables, and in terms of arity of constraints) 
in spite of the exponential complexity of the algorithms. These implementations 
generate CHR rules that deal with finite domain variables using an ECL'PS'^ library. 
Constraint Handling Rules (CHR) of ( Friihwirth, 1995| ) is a declarative language 



that allows one to write guarded rules for rewriting constraints. These rules are re- 
peatedly applied until a fixpoint is reached. The rule applications have a precedence 
over the usual resolution step of logic programming. 

A CHR program is a finite set of CHR rules. These rules are basically of two types 
(there is a third type of rules which is a combination of the first two types): sim- 
plification rules and propagation rules. When all guards are satisfied, a simplifica- 
tion rule replaces constraints by simpler ones while preserving logical equivalence, 
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whereas a propagation rule adds logically redundant constraints. More precisely, 
these rules have the following form: 

simplification Hi, . . . ,Hi <=> Gi, . . . ,Gj \ Bi, . . . ,Bk 
propagation Hi,..., Hi Gi,...,Gj \ Bi, . . . , Bk 

where 

• j > 0, j > 0, fc > 0, 

• the multi-head Hi, . . . ,Hi is a non-empty sequence of CHR constraints, 

• the guard Gi, . . . ,Gj is a sequence of built-in constraints, 

• the body Bi, . . . , Bk is a. sequence of built-in and CHR constraints. 

Our equality rules and membership rules can be modelled by means of propaga- 
tion rules. To illustrate this point consider some constraint cons on three variables, 
A, B, C, each with the domain {0, 1, 2}. 

The Rules Generation algorithm generates rules such as iA,C) = (0,1) — » 
B ^2. This rule is translated into a CHR rule of the form: cons(0,B, 1) ==> B##2. 
Now, when a constraint in the program query unifies with cons(0,B,l), this rule 
is fired and the value 2 is removed from the domain of the variable B. 

In turn, the MEMBERSHIP Rules Generation algorithm generates rules such 
as (A, G) e ({0}, {1, 2}) ^ S ^ 2. This rule is translated into the CHR rule 

cons(0,B,C) ==>in(C, [1,2] ) I B##2 

where the in predicate is defined by 

in(X,L):- dom(X,D), subset(D,L). 

So in(X,L) holds if the current domain of the variable X (yielded by the built-in 
dom of ECL*PS'=) is included in the list L. 

Now, when a constraint unifies with cons(0,B,C) and the current domain of the 
variable C is included in [1,2], the value 2 is removed from the domain of B. So for 
both types of rules we achieve the desired effect. 

In the examples below we combine the rules with the same premise into one rule 
in an obvious way and present these rules in the CHR syntax. 



10.2 Generating the rules 

We begin by discussing the generation of equality rules and membership rules for 
some selected domains. The times given refer to an implementation ran on a Silicon 
Graphics 02 with 64 Mbytes of memory and a 180 MHZ processor. 



Boolean constraints As the first example consider the Boolean constraints, for 
example the conjunction constraint and(X,Y,Z) of Example 5.2. The Equality 
Rules Generation algorithm generated in 0.02 seconds the following six equality 
rules: 



and(l,l,X) ==> X##0. 
and(X,0,Y) ==> Y##l. 
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and(0,X,Y) ==> Y##l. 
and(X,Y,l) ==> X##0,Y##0. 
and(l,X,0) ==> X##l. 
and (X, 1,0) ==> X##l. 

Because the domains are here binary we can replace the conclusions of the form 
U ## by U = 1 and U ## 1 by U = 0. These rules are somewhat different than 
the well-known rules that can be found e.g. in ( Friihwirth, 1998| )[page 113], where 
instead of the rules 

and(l,l,X) ==> X##0. 
and(l,X,0) ==> X##l. 
and(X,l,0) ==> X##l. 

the rules 

and(l,X,Y) ==> X = Y. 
and(X,l,Y) ==> X = Y. 
and(X,Y,Z), X = Y ==> Y = Z . 

appear. We shall discuss this matter in Section |ll|. 

In this case, by virtue of Theorem ■iA, the notions of rule and arc consistency 
coincide, so the above six equality rules characterize the arc consistency of the and 
constraint. Our implementations of the Equality Rules Generation and the 
Membership Rules Generation algorithms yield here the same rules. 



Three valued logic Next, consider the three valued logic of ( Kleene, 1952[ )[page 334] 
that consists of three values, t (true), f (false) and u (unknown). We only consider 
here the crucial equivalence relation = defined by the truth table 







f 


u 


t 


t 


f 


u 


f 


f 


t 


u 


u 


u 


u 


u 



that determines a ternary constraint with nine triples. We obtain for it 20 equality 
rules and 26 membership rules. Typical examples are 

equiv(X,Y,f) ==> X##u,Y##u. 

and 

equiv(t,X,Y) ==> in(Y,[f, u] ) I X##t . 



Six valued logic In ( Van Hcntcnryck et a/., 1992| ) the constraint logic programming 
language CHIP is used for the automatic test-pattern generation (ATPG) for the 
digital circuits. To this end the authors define a specific six valued logic and provide 
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some rules (expressed in the form of so-called demons) to carry out the constraint 
propagation. The andS constraint in question is defined by means of the following 
table: 



and6 1 d dnot e enot 












1 
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d 
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dnot 




dnot 
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enot 





enot 



— — 
d dnot c enot 
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dnot 

d e 

dnot enot 



The Equality Rules Generation algorithm generated 41 equality rules in 
0.15 seconds, while the Membership Rules Generation algorithm generated 
155 membership rules in 14.35 seconds. The generated rules enforce, respectively, 
rule consistency and arc consistency, while it is not clear what notion of local 
consistency is enforced by the (valid) rules of ( Van Hcntcnryck et ai, 1992| )[page 
133] because some of the latter ones allow equalities between the variables in the 
premise. This makes the comparison in terms of strength of the entailed notion of 
local consistency difficult. It is clear that our approach is more systematic and fully 



automatic. (In fact, we found two typo's in the rules of (Van Hentenryck et 



1992| )[page 133].) 



Propagating signs As a next example consider the rules for propagating signs in 
arithmetic expressions, see, e.g., (Davis, 1987) [page 303]. We limit ourselves to the 
case of multiplication. Consider the following table: 



X neg zero pos unk 



neg pos zero neg unk 

zero zero zero zero zero 

pos neg zero pos unk 

unk unk zero unk unk 



This table determines a ternary constraint msign that consists of 16 triples, for 
instance (neg, neg, pos) that denotes the fact that the multiplication of two negative 
numbers yields a positive number. The value "unk" stands for "unknown" . The 
Equality Rules Generation algorithm generated in 0.08 seconds 34 equality 
rules. A typical example is msign (X, zero, Y) ==> Y##pos , Y##neg, Y##unk. 

In turn, the MEMBERSHIP Rules Generation algorithm generated in 0.6 sec- 
onds 54 membership rules. A typical example is 

msign (X, unk, Y) ==> in (Y, [neg, pos, zero]) I X##pos,X##neg 
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that corresponds to the following two membership rules for the constraint msign (X , Z , Y) : 
(Z, Y) G ({unk}, {neg, pos, zero}) X ^ pos 

and 

{Z, Y) G ({unk}, {neg, pos, zero}) X ^ neg. 



Waltz' language for describing polyhedral scenes Waltz' language consists of four 



constraints. One of them, the fork junction was already mentioned in Example 9.1 



The Equality Rules Generation algorithm generated for it 12 equality rules 
and the Membership Rules Generation algorithm 24 membership rules. 
Another constraint, the so-called T junction, is defined by the following table: 



X I y 



In this case the Equality Rules Generation algorithm and the Membership 
Rules Generation algorithm both generate the same output that consists of just 
one rule: 

t(X,Y,Z) ==> X##'l' ,X##'-' ,X##'+' ,Y##'r' ,¥##'-',¥##'+' . 

So this rule characterizes both rule consistency and arc consistency for the CSP's 
based on the T junction. 

For the other two constraints, the L junction and the arrow junction, the gener- 
ation of the equality rules and membership rules is equally straightforward. 



10.3 Using the rules 

Next, we show by means of some examples how the generated rules can be used to 
reduce or to solve specific queries. Also, we show how using compound constraints 
we can achieve local consistency notions that are stronger than arc consistency for 
constraints of arbitrary arity. 

Waltz' language for describing polyhedral scenes The following predicate describes 
the impossible scene given in Figure |l] and taken from ([Winston, 1992|) [page 262]: 



imp(AF,AI,AB,IJ,IH, JH,GH,GC,GE,EF,ED,CD,CB) :- 

S1=[AF,AI,AB,IJ,IH, JH,GH,GC,GE,EF,ED,CD,CB] , 
S2=[FA,IA,BA, JI,HI,HJ,HG,CG,EG,FE,DE,DC,BC] , 
append(Sl,S2,S) , S :: [+,-,l,r], 




arrow ( AF, AB, AI ) , 1(BC,BA), arrow (CB, CD, CG) , 1(DE,DC), 
arrow (ED, EG, EF) , 1(FA,FE), f ork(GH,GC,GE) , arrow(HG,HI ,HJ) , 
fork(IA,IJ,IH) , 1(JH,JI), 

line(AF,FA) , line(AB,BA), line (AI , lA) , line(IJ,JI), 
liiie(IH,HI) , line(JH,HJ), liiie(GH,HG) , line(FE,EF), 
line(GE,EG), line(GC,CG), line(DC,CD), line(ED,DE), 
line(BC,CB) . 

where the supplementary constraint line is defined by the following table: 



|x|y| 

+ + 

1 r 
r 1 



Here and elsewhere we use the ECL^PS*^ built-in : : to declare variable domains. 

When using the equality rules obtained by the Equality Rules Generation 
algorithm and associated with the fork, arrow, t, 1, and line constraints, the 
query 

imp(AF,AI,AB,IJ,IH,JH,GH,GC,GE,EF,ED,CD,CB) 

reduces in 0.009 seconds the variable domains to AF € [+,-, 1] , AI € [+,-], AB 
e [+,-,r],IJ G [+,-,l,r],IH G [+,-,l,r],JH € [+,-,l,r], 
GH € [+,-,l,r],GC S [+,-,l,r],GE € [+,-,l,r],EF G [+,-], 
ED e [+,-,l],CD e [+,-,r] , and CB G [+,-,1]. 

But some constraints remain unsolved, so we need to add a labeling mechanism 
to prove the inconsistency of the problem. On the other hand, when using the mem- 
bership rules, the inconsistency is detected without any labeling in 0.06 seconds. 
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In the well-known example of the cube given in Figure 12.15 of ( [Winston, 199^ ) [page 
260] the membership rules are also more powerful than the equality rules and both 
sets of rules reduce the problem but in both cases labeling is needed to produce all 
four solutions. 

Comparing the constraint solver based on the membership rules to the constraint 
solver based on the equality rules is not easy: although propagation is more efficient 
with the membership rules, the solver based on the equality rules can sometimes 
be faster depending on the structure of the problem and on whether the labeling is 
needed. 

We also compared the solvers generated by the implementations of the Equality 
Rules Generation and Membership Rules Generation algorithms to the 
approach described in (By, 1997) and based on meta-programming in Prolog. We 
ran the same examples and drew the following conclusions. For small examples 
our solvers were less efficient than the ones of By (with factors varying from 2 to 
10). However, for more complex examples, our solvers became significantly more 
efficient, with factors varying from 10 to 500. This can be attributed to the increased 
role of the constraint propagation that reduces backtracking and that is absent in 
By's approach. 



Temporal reasoning In (Allen, 1983 )'s approach to temporal reasoning the entities 



are intervals and the relations are temporal binary relations between them. ( Allen 



1983) found that there are 13 possible temporal relations between a pair of events, 
namely before , during, overlaps, meets, starts, finishes, the symmetric 
relations of these six relations and equal. We denote these 13 relations respectively 
by b,d, o ,m, s ,f ,b- ,d- , o- ,in- , S-, f - ,e and their set by TEMP. 

Consider now three events. A, B and C and suppose that we know the temporal 
relations between the pairs A and B, and B and C. The question is what is the 



temporal relation between A and C. To answer it (Allen, 1983) provided a 13 x 13 
table. This table determines a ternary constraint between a triple of events. A, B 
and C that we denote by alien. For example, 

(overlaps, before, before) G alien 

since A overlaps B and B is before C implies that A is before C. 

Using this table, the Equality Rules Generation algorithm produced for the 
constraint alien 498 equality rules in 31.16 seconds. In contrast, we were unable 
to generate all membership rules in less than 24 hours. This shows the limitations 
of our implementation. We tried the generated set of equality rules to solve the 



following problem from (Allen, 1983): "John was not in the room when I touched 
the switch to turn on the light." . We have here three events: S, the time of touching 
the switch; L, the time the light was on; and J, the time that John was in the room. 
Further, we have two relations: Rl between L and S, and R2 between S and J. This 
problem is translated into the CSP (alien ; Rl e [o— ,m— ],R2 £ [b, m, b— , m— ], R3 G 
TEMP), where alien is the above constraint on the variables Rl, R2, R3. 
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To infer the relation R3 between L and J we can use the foUowing query|^ 

Rl : : [o- ,m-] , 

R2 : : [b,m,b-,m-] , 

R3 : : [b,d,o,m,s,f ,b-,d-,o-,m-,s-,f-,e] , 
alien (R1,R2,R3) , 
labeling ( [R1,R2,R3] ) . 

We then obtain the foUowing solutions in 0.06 seconds: 

(R1,R2,R3) e {(m-,b,b) , (m-,b,d-), (m-,b,f-), (m-,b,m) , (m-,b,o), 
(m-,b-,b-), (m-,m,e) , (m- , m , s ) , (m- ,m, s-) , (m-,m-,b-), (o-,b,b), 
(o-,b,d-), (o-,b,f-), (o-,b,m), (o-,b,o), (o-,b-,b-), (o-,m,d-), 
(o-,m,f-), (o-,m,o), (o-,m-,b-)}. 



To carry on (as in (Allen, 1983)), we now complete the problem with: "But 
John was in the room later while the light went out.". This is translated into: "L 
overlaps, starts, or is during J", i.e., R3 G [o,s,d]. 

We now run the following query: 

Rl : : [o- ,m-] , 
R2 : : [b,m,b-,m-] , 
R3 : : [o,s,d] , 
allen(Rl,R2,R3) , 
labeling ( [R1,R2,R3] ) . 

and obtain four solutions in 0.04 seconds: (R1,R2,R3) G {(m-,b,o), (m-,m,s), 
(o-,b,o) , (o-,m,o)}. 

Three valued logic Next, consider the and3 constraint in the three valued logic of 



(Kleene, 1952) [page 334] represented by the truth table 
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Typical examples of the 16 generated equality rules and of the 18 generated 
membership rules are: 

and3(u,u,X) ==> X##0 . 

and 

and3(X,Y,Z) ==> in(X,[0, u] ) I Z##l. 

Consider now the query: 



^ Since no variable is instantiated, we need to perform labeling to effectively apply the rules. 
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[X,Y,Z,T,U] : : [0,l,u] , and3(X,Y,Z) , and3(T,U,Z), Z##0,Y##u, X##u. 

Using the membership rules we reach the answer in the form of a complete as- 
signment to all the variables, namely X = 1,Y = 1,Z = 1,T = 1,U = 1, whereas 
using only the equality rules we do not obtain any reduction and labeling is needed 
to produce the same answer. 



Full adder It is often the case that dealing with a compound constraint directly 
yields a stronger notion of local consistency than when dealing with each of its 
constituents separately. A prototypical example is the well-known alldif f erent 
constraint on n variables that can be decomposed into binary disequality 

constraints. Then arc consistency enforced on alldif f erent is stronger than arc 
consistency enforced on each of the disequality constraints separately. 

The same phenomenon arises for rule consistency. We illustrate it by means 



of the already discussed in Example 6.3 full adder circuit. It can be defined by 



the following constraint logic program (see, e.g., ( Friihwirth, 1998 )) that uses the 
Boolean constraints and, xor and or: 

add(Il,I2,I3,01,02) :- 

[I1,I2,I3,01,02,A1,A2,X1] : : 0. .1, 

xor(Il,I2,Xl) , 

and(Il,I2,Al) , 

xor(Xl,I3,02) , 

and(I3,Xl,A2) , 

or(Al,A2,01) . 

The query add(Il,I2,I3,01,02) followed by a labeling mechanism generates the 
explicit definition (truth table) of the full_adder constraint with eight entries such 
as 

full_adder( 1,0, 1,1,0) . 

We can now generate the equality rules for the compound constraint (here the 
full_adder constraint) that is defined by means of some basic constraints (here 
the and, or and xor constraints). These rules refer to the compound constraint 
and allow us to reason about it directly instead of by using the rules that deal with 
the basic constraints. 

In the case of the f ull_adder constraint the Equality Rules Generation 
algorithm generated 52 equality rules in 0.27 seconds. The constraint propagation 
carried out by means of these equality rules is more powerful than the one carried 
out by means of the equality rules generated for the and, or and xor constraints. 
For example, the query [X,Y,Z] : : [0,1] , full_adder(l,X,Y,Z,0) reduces Z to 
1 whereas the query [X,Y,Z] : : [0,1] , add(l ,X,Y,Z,0) does not reduce Z at all. 

So rule consistency for a compound constraint defined by means of the basic 
constraints is indeed in general stronger than the rule consistency for the basic 
constraints treated separately. In fact, in the above case the equality rules for the 



full_adder constraint yield the relational (l,5)-consistency notion of (Dcchter 
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van Beek, 1997), whereas by virtue of Theorem 8.4, the equahty rules for the and, 



or and xor constraints yield a weaker notion of arc consistency. 



11 Related Work 

11.1 Relation Between Constraint Programming and Rule-Based 

Programming 

In a number of papers a link was made between constraint programming and rule- 



based programming. To start with, in (Montanari & Rossi, 1991) a general study 



of constraint propagation was undertaken by defining the notion of a relaxation 
rule and by proposing a general relaxation algorithm that implements constraint 
propagation by means of a repeated application of the relaxation rules. However, 
this abstract view of constraint programming cannot be realized in a simple way 
since the application of a relaxation rule is a complex process. 



In (Apt, 199S) we showed how constraint programming can be couched in proof 
theoretic terms by viewing the programming process as the task of proving the 
original CSP. In the proposed framework two types of rules were proposed: deter- 
ministic ones and the splitting ones. Further, the deterministic rules were either 
concerned with domain reduction or constraint reduction. In the former case the 
rules were called domain reduction rules and in the latter case constraint reduction 
rules. Such rules are high-level abstractions and on the implementation level they 
can involve complex computations. 

It is useful to see that the rule-based approach to constraint programming pro- 
posed in this paper is an instance of this proof theoretic view of constraint pro- 
gramming. Namely, the equality rules and the membership rules are examples of 
the domain reduction rules while labeling, the formal treatment of which is omitted 
here, is an example of a splitting rule. 

The important gain is that the implementation of the considered here equality 
rules and membership rules boils down to a straightforward translation of them 
into the CHR syntax. This leads to an implementation of this approach to constraint 
programming by means of constraint logic programming. The important limitation 
is that this approach applies only to the CSP's built out of predefined, explicitly 
given finite constraints. 



A similar approach to constraint programming to that of (Apt, 1998) was pro- 



posed in (Castro, 199S). In his approach the proof rules are represented as rewrite 
rules in the already mentioned in the introduction programming language ELAN. 
The rules use a richer syntax than here considered by referring to arbitrary con- 
straints and to expressions of the form x G D, where D is the current domain of the 
variable x. In particular no constraint specific rules were considered. Instead, the 
emphasis was on showing how the general techniques of constraint programming, 
in particular various search strategies, can be expressed in the form of rules. 
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11.2 Generation of Rules 



Let us turn now to an overview of the recent work on rules generation. Building 



upon the work presented in (Apt & Monfroy, 1999) two articles appeared in which 
algorithms were presented that aim at improving the expressivity of the rules and 
at a more economic representation. 



In (Ringeissen & Monfroy, 2000) rules similar to equality rules were considered. 
The most significant improvement is the use of parameters (i.e., unspecified con- 
stants) that leads to a decrease in the number of generated rules. Parameters are 
also a means for deducing equalities of variables in the right-hand side of rules. For 
instance, consider the following two rules with parameters ai and 02 for a constraint 



C taken from ( [Ringeissen fc Monfroy, 2000[) 



xi = ai 



X2 



Oi A = 1 



(1) 



a;2 = a2 xi = 02 A 2:3 = 1 (2) 

Rule (^ means that whatever the value of xi is, X2 is equal to xi and 2:3 is equal to 
1. From rules (|l|) and (||) an equality between the variables on the right-hand side 
of the rule can be deduced (note that the resulting rule can always be applied): 

^ a;2 a;i A a;3 = 1 



To generate such rules with parameters ( [Ringeissen fc Monfroy, 200C ) combine 



unification in finite algebra with a rule generation algorithm. The size of the gener- 
ated set of rules significantly depends on an ordering on variables. It still needs to be 
clarified what is the counterpart of the notion of a minimal rule in this framework 
and whether the generated rules with parameters enforce rule-consistency. 



Sets of rules generated in (Abdennadher & Rigotti, 2000) are even more compact 



and more expressive: multiple occurrences of variables and conjunction of con- 
straints with shared variables are allowed in the left-hand side of rules. Moreover, 
the user has the possibility to specify the admissible syntactic form of the rules: 
more specifically, right-hand sides of rules can consists of more complex constraints 
than (dis)-equality constraints. Here are two examples of rules (taken from two 
different sets of rules of ( Abdennadher fc Rigotti, 2000| )): 



and{x, x,z) ^ X = z. (3) 
and{x,y,z),neg{x,y)-^z — Q. (4) 

In rule ^ equality between variables is deduced using a double occurrence of the 
variable x in the head, and rule (^) defines interaction between two constraints, 
and and neg. 



(Abdennadher & Rigotti, 2000) also investigated what form of local consistency is 
enforced by the rules generated by their algorithms. When using the given constraint 
together with equality (between a variable and a value) on the left-hand side, and 
only disequality (between a variable and a value) on the right-hand side of the 
rules, the generated rules enforce rule consistency. However, in general, the enforced 
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local consistency is stronger than rule consistency. In particular, it is plausible that 
membership rule consistency (i.e., arc consistency) is enforced when disequality 
constraints are allowed both on the right- and left-hand sides of rules. 



In (Apt, 2000) two sets of rules for Boolean constraints were compared. One of 



them is the one presented in Subsection 10.2 and the other the already mentioned 



set of rules from (Friihwirth, 1998) [page 113] (with one difference irrelevant for the 
subsequent discussion). While both sets of rules enforce the arc consistency, it turns 
out that they are not equivalent. In fact, if a Boolean CSP with non-empty domains 
is closed under the rules from the second set, then it is closed under the rules from 
the first set. The converse does not hold since the CSP {x A y — z ; x & {1},?/ 6 
{0, 1}, z £ {0, 1}) is closed under the first set of rules but not under the second one, 
since it is not closed under the rule and(l,X,Y) ==> X = Y. 



11.3 Local Consistency Notions 

In this paper we considered local consistency by focusing on individual, arbitrary, 
non-binary, constraints. In the literature algorithms for achieving such notions of 
local consistency usually concentrated on achieving arc consistency of non-binary 



CSP's. Among them the algorithm CN of (Mackworth, 1977b), and GAC4 of (Mohr 



|fc Masini, 1988 ) are respectively based on ideas similar to the AC-3 and AC-4 



algorithms for binary constraints. However, in practice CN is usable only for ternary 
constraints and small domains, and has a large worst-case time complexity. On the 
other hand, the large space complexity of GAC4 makes it usable only for constraints 
of small size. 



The GAC- schema of ( Bessiere &: Regin, 1997 ) was designed to enforce arc consis 



tency on non-binary constraints while keeping a reasonable time and space com- 
plexity. It is based on an AC-7 like schema and allows constraints to be given 
explicitly, either in a positive way as in our case, or in a negative way, i.e., in the 
form of "forbidden" tuples, or implicitly in the form of predicates. In order to make 
our framework as general as the GAC- schema, we could think of generating the al- 
lowed tuples by testing all possible tuples. However, this would almost always be 
impractical because of space considerations. 



In (Bessiere, 1999) the following opinion was voiced on non-binary constraints: 
"Perhaps we should accept the idea that the constraint solving tool of the next 
years will apply different levels of local consistency on different constraints at each 
node of the search tree, Our framework is amenable to such a view since we 
can generate rules for enforcing various types of local consistencies. For example, 
we can generate some equality rules for some constraints and some membership 
rules for other constraints, and then apply only the resulting set of rules. 

Another consideration is that our constraint solving process consists of two sep- 
arate phases: first a generation of the rules (a sort of compilation of the truth table 
of a constraint) which is done once and for all, followed by the application of the 
rules. Thus, the set of generated rules can be modified during the application phase, 
for example by combining some rules (for a more efficient application of the rules), 
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by removing some rules or by strengthening some conditions to weaken the domain 
reduction process. 

While preparing this revised version of the paper we noted that a similar notion 
to our rule consistency notion was introduced in the context of the theory of fuzzy 



sets, see (Pedrycz fc Gomide, 1998D [pages 252-261]. The notion there considered 



deals with rules of the form "if x \s A, then y is B" , where A and B are fuzzy sets. 
In spite of the same name used (namely, rule consistency) , the uses of both notions 
are different. In our case, we employ it to reduce a specific CSP to a smaller one that 
is rule consistent but can be inconsistent, while in the case of the fuzzy set theory 
the corresponding notion is used to detect conditions for "potential inconsistency" 
that arises when the rules express contradictory knowledge. 

Finally, let us mention that rule generation appears in other areas of computer 
science. Typical examples are: programs for machine learning that construct a model 



of the knowledge using decision trees and production rules (see, e.g., ( Quinlan, 



1993 )) inductive logic programming, a logic-based approach to machine learning 



where logic programming rules are inferred from positive and negative examples 



and a background knowledge (see, e.g., (Muggleton fc de Raedt, 1994)), and data 



mining that aims at extracting high-level representations in the form of patterns 



and models from data (see, e.g., ( Agrawal et al, 1996 ) where so-called association 
rules are generated.). 



12 Conclusions 

The aim of this paper was to provide a framework in which constraint program- 
ming can be entirely reduced to rule-based programming. It involved constraint 
satisfaction problems built out of explicitly given constraints. In the case the latter 
constraints are defined over small finite domains these CSP's can be often solved 
by means of automatically generated constraint propagation algorithms. 

We argued that such CSP's often arise in practice and consequently the methods 
here developed can be of practical use. We believe that the approach of this paper 
could be applied to a study of various decision problems concerning specific multi- 
valued logics and this in turn could be used for an analysis of digital circuits (see. 



e.g., ( Muth, 1976 ) where a nine valued logic is used). Other applications could in- 
volve non-linear constraints over small finite domains and the analysis of polyhedral 
scenes in presence of shadows (see ( Waltz, 1975| )). 



The introduced notion of rule consistency is weaker than arc consistency and can 
be in some circumstances the more appropriate one to use. For example, for the 
case of temporal reasoning considered in the last section we easily generated all 
498 equality rules that enforce rule consistency whereas 24 hours turned out not be 
enough to generate the membership rules that enforce arc consistency. (For a more 
precise summary of the tests carried out see the appendix.) 

In this paper, we focused on systematic and automated aspects of rule-based 
constraint solvers. At present stage it is difficult to compare the performance of 
our method (based on rule generation and subsequent rule application) with other 
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methods based on classical constraint propagation algorithms. The reason is that 
our approach is currently implemented by means of CHR rules that are applied 
on top of Prolog while the built-in constraint propagation algorithms are usually 
implemented at a lower level. 

Finally, the notions of rule consistency and membership rule consistency could be 
parametrized by the desired maximal number of variables used in the rule premises. 
Such parametrized versions of these notions could be useful when dealing with con- 
straints involving a large number of variables. Both the Equality Rules Gener- 
ation algorithm and the Membership Rules Generation algorithm and their 
implementations can be trivially adapted to such parametrized notions. 

The approach proposed in this paper could be easily integrated into constraint 
logic programming systems such as ECL'PS'^. This could be done by providing an 
automatic constraint propagation by means of the equality rules or the membership 
rules for flagged predicates that are defined by a list of ground facts, much in the 
same way as now constraint propagation for linear constraints over finite systems 
is automatically provided. 
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Appendix 

Table ^ illustrates the generation of the equality rules and the membership rules for 
various natural constraints. The first column gives the name of the constraint, the 
second its arity, the third the cardinality of the domains of variables, and the fourth 
the cardinality of the constraint. The subsequent two columns show the outcome 
of the Equality Rule Generation algorithm: first the number of equality rules 
generated, and then the computation time in seconds; the last two columns provide 
this information for the membership rules. 

The fork, t, msign, alien, and fulladder constraints represent the previously de- 
scribed constraints. The noti and andi constraints, where i £ {3, 4, 6, 8, 9}, represent 
the usual not and and operators for multi- valued logics. The blOm constraint is the 
multiplication of digits from to 9, i.e., blOm{X, Y, C, Z) stands for the constraint 
X*y = Z -\- IQ* C defined over the intervals [0..9]. The "— " symbol means that 
we were unable to generate the rules in less than 24 hours. 

The constraints for multi-valued logics are presented in a way that shows the 
impact of the domain size and of the cardinality of the constraint in case of the 
same arity and a similar structure. 

In spite of its exponential running time, the Equality Rules Generation 
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Constraint 


Aritv 


Domain 


Cctrdinality 


Equality 


Gen. 


Member. 


Gen. 






Size 




Rules 


fin s ^ 


Rules 


fin s 1 


fork 


3 


4 


5 


12 


0.05 


24 


0.65 


t 


3 


4 


4 


1 


0.02 


1 


0.07 


not 


2 


2 


2 


4 


0.01 


4 


0.19 


nots 


2 


3 


3 


6 


0.02 


6 


0.29 


noti 


2 


4 


4 


8 


0.02 


8 


0.5 


note 


2 


6 


6 


12 


0.03 


12 


0.14 


nots 


9 


a 

o 


fi 

o 


1 fi 


u.uo 


1 R 


n fi7 


notg 


2 


9 


9 


18 


0.07 


18 


1.57 


and 


3 


2 


4 


6 


0.02 


6 


0.08 


ands 


3 


3 


9 


16 


0.04 


18 


0.13 


andi 


3 


4 


16 


26 


0.08 


43 


0.6 


ande 


3 


6 


24 


41 


0.15 


155 


14.35 


ands 


3 


8 


64 


96 


0.57 


622 


351.16 


andg 


3 


9 


81 


134 


1.07 


1294 


1777 


msign 


3 


4 


16 


34 


0.08 


54 


0.6 


fulladder 


5 


2 


8 


52 


0.29 


52 


0.38 


610m 


4 


10 


100 


362 


14.83 






alien 


3 


13 


409 


498 


31.16 







Table 1. Generation of equality and membership rules 

algorithm is stiU usable. On the other hand, the Membership Rules Generation 
algorithm is much more costly for larger problems. Sometimes it also generates too 
many rules for medium size problems (such as andg) and thus becomes unusable. 
In general, it is difficult to decide which notion of local consistency should be used 



to solve a given CSP. In particular (Sabin & Freudcr, 1994) showed that in the case 



of CSP's consisting of binary constraints maintaining full arc consistency during the 
backtracking search can be often more efhcient than a more limited of constraint 
propagation embodied in the so-called forward checking. However, empirical results 
for CSP's involving non-binary are missing and it is quite conceivable that for 
such CSP's imposing full arc consistency during the backtracking search can be 
too costly. For these CSP's a weaker form of constraint propagation, such as rule 
consistency, could be an alternative. 



